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A systematic connection between QCD and nuclear few- and many-body properties in the form 
of the Effective Field Theory "without pions" is applied to A < 6 nuclei to determine its range of 
applicability. We present results at next-to-leading order for the Tjon correlation and for a corre- 
lation between the singlet S-wave 3 He-neutron scattering length and the triton binding energy. In 
the A = 6 sector we performed leading order calculations for the binding energy and the charge 
and matter radii of the halo nucleus 6 He. Also at leading order, the doublet S-wave 4-He-neutron 
phase shifts are compared with R-matrix data. These analysis provide evidence for a sufficiently 
fast convergence of the effective field theory, in particular, our results in A < 4 predict an expan- 
sion parameter of about i, and they converge to data within the predicted uncertainty band at this 
order. A properly adjusted three-body contact force which we include together with the Coulomb 
interaction in all calculations is found to correctly renormalize the pion-less theory at leading- and 
next-to-leading order, i.e., the power counting does not require four-body forces at the respective 
order. 
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1. Introduction 

In principle, effective field theories (EFT) provide a framework in which properties of nuclei 
can be derived from the underlying theory of QCD. In contrast to the elegance and success of Chiral 
Perturbation theory in the mesonic sector, its extention, by including baryonic degrees of freedom, 
called Chiral Effective Field Theory, suffers practical and systematic problems if it is used to ex- 
pand multi-baryon amplitudes [jl|]. There are no such inconsistencies in the power counting of the 
Effective Field Theory "without pions" (EFT^, see e.g. Furthermore, this EFT has a simple 
vertex structure with nucleon-nucleon (NN) contact terms only, which makes calculations in larger 
nuclei more practical. 

EFT^ systematically expands observables in powers of the typical low-momentum scale p typ of 
the process considered, measured in units of the scale A/, at which the pion can be resolved as a 
dynamical degree of freedom, set by the pion mass m n . While this expansion has been shown to 
converge sufficiently fast for many reactions in the two- and three-nucleon sector with and without 
electro-weak currents (see e.g. its usefulness for systems with A > 4 has been under question. 
The question whether or not an observable is apt to be expanded within EFT^ is relatively easy to 
answer for two nucleons where p typ is closely related to, e.g., the relative momentum in neutron- 
proton scattering, or the binding momentum of the deuteron. In contrast, the 4 He binding energy 
B( 4 He) = 28.5 MeV, and its typical size of about 1.6 fm is comparable with the approximate range 
of the one-pion exchange. A systematic expansion in the four-nucleon system could therefore con- 
verge very slowly, at best, if the naive estimate for p typ turns out to be appropriate. The current 
understanding of effective field theories does not provide rigorous arguments at which values, e.g., 
of A or of the density inside the nucleus, EFT^ fails. Explicit calculations are required! 
To demonstrate the usefulness of EFT^ applied to the 4-nucleon bound- and scattering system, we 
analyzed the Tjon correlation and a correlation between the singlet S-wave 3 He-neutron scattering 
length ao( 3 He-n) and the triton binding energy B(t). The results are taken from the more detailed 
analysis in [Q] (an extensive list of references can be found there). With the systematic and predic- 
tive power thus demonstrated, the status of our investigation in an A = 5 scattering system and the 
halo system 6 He is shown. Specifically, we present leading order (LO) results for the doublet S- 
wave phase shift for 4 He-neutron scattering, in addition to binding energy, charge- and matter radii 
results for 6 He. Before that, we present the specific form of the next-to-leading order (NLO) EFT^ 
NN interaction and the few -body technique of the Refined Resonating Group Model (RRGM ^) 
used to solve the many-body Schrodinger equation. 



2. Pion-less Theory of the Nuclear Interaction 

The interaction part of the Lagrangean of EFT^ at NLO (see e.g. [§]) consists of nine four- 
nucleon contact terms and one six-nucleon vertex. A spin dependend term, one spin independent 
one, and the six-nucleon term are momentum independent and comprise the LO part of the potential 
which is to be iterated to all orders because of the unnaturally large scattering lengths a s j in the 
NN system. At NLO vanishing P-wave amplitudes result in dependencies amongst the seven new 
coupling strengths and reduce their number to only three independent NLO parameters. 
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We solve the Schrodinger equation in coordinate space with the following operator basis: 
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(2.1) 



By changing the renormalization scheme in two distinct ways, we test the consistency of the power 
counting and get an estimate of higher order effects, i.e., the convergence rate of the theory. The 
two methods are: 

1. A variation of the parameter A of the Gaussian regulator in the interval [400; 700] MeV. 

2. A variation of the data used to fit the low-energy constants (LEC) A,-: We did not implement 
the dependencies between the A, 's but fitted nine coupling strengths independently to allow 
for non-zero P-wave phase shifts up to 0.1% of the Nijmegen values. 

Both represent a modification of the unresolved short distance structure of the nucleons and should 
affect low-energy observables only within the uncertainty range of the considered order. Nine 
NLO potentials, corresponding to nine different sets of A,'s, were fitted to the deuteron binding 
energy B(d) and to the neutron-proton (np) phase shifts d(^' 3 'So t i) , £1 ( 3 5i — 3 D\ ) below 0.3 MeV 
from |^]. Each LEC set was determined either with a different cutoff A or different values for 
the P-wave phase shifts. The three-body input to fix At,nf was B(t), and the LO parameters were 
adjusted to the singlet and triplet np scattering lengths. 

Corrections to the non-relativistic kinetic energy term are of higher order. Hence, amplitudes are 
calculated by first iterating the LO part of the potential to all orders by solving, e.g. , the Schrodinger 
equation, and a consecutive perturbative treatment of the NLO operators in distorted wave Born 
approximation. However, we iterate also the NLO vertices and conclude from the results that 
higher order diagrams which are included in such an approach yield only negligible contributions. 

3. Resonating Group Model 

To reach the goals stated in the introduction, we need a method to solve the few-body bound- 
and scattering problem. We employ a refined version of the Resonating Group Model which yields 
a variational solution to the Schrodinger equation in coordinate representation (see e.g. The 
variational space for a N-body bound state is spanned by a superposition of antisymmetric 
states of the form: 



, dJJ 



k=\ 



(3.1) 



The index labels different (orbital) spin angular momentum coupling schemes and thus resem- 
bles the idea of the resonating groups, i.e., a hierarchy of all possible groupings of neutrons and 
protons within a nucleus which allows a simplification of the wave function ansatz by discarding 
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those states corresponding to groupings that are relatively short lived. This simplification, irrel- 
evant thanks to modern computer technology for the three- and four-body observables for which 
we included all possible groupings, becomes an efficient tool that makes A > 5 body calculations 
feasible. For instance, a single fragmentation with a 4 He core and two attached neutrons was found 
to yield more than 95% of the assumed converged observables calculated here for 6 He. 
The dimension of the model space is further reduced by including only orbital angular momenta 
hi < 2 on the Jacobi coordinates p. In addition, certain couplings even with 1 = 1, e.g., for the 
coordinate between the 4 He core and the center of mass of the two skin neutrons and the relative 
coordinate between those neutrons, an explicit calculation showed the contribution to B( 6 He) to be 
negligible. 

Finally, the number and magnitudes of the Gaussian width parameters y c / had to be tailored to the 
various potentials to allow for enough flexibility of the basis in the range where the potential is non- 
zero. This interval varies with the cutoff A. The numerical instabilities due to the non-orthogonality 
of the basis were monitored to ensure also a numerically independent set of basis vectors. 

4. Results 

4.1 Four nucleons: bound state 

EFT^ should appropriately describe the 4 He bound state to make its application to larger sys- 
tems, where this state plays an important role as a target or core, a worthwile endeavor. We analyzed 
the dependence of B( 4 He) on B(t) using nine NLO potentials with an unconstrained Aj^p param- 
eter, i.e., the two-body sector is renormalized properly but for A > 2 one expects deviations from 
data analogous to older model calculations with NN-forces only. With an EFT, this dependence, 
known as Tjon-line, becomes a band whose width is a measure for the accuracy of the expansion. 
A change in the renormalization scheme, in our case a change in the cutoff A and the P-wave 
phase shifts, should yield results for the binding energy tuple within an uncertainty range set by 
(p tyv /m n ) n around the datum in an N"LO calculation. 

In the left panel of Fig. [I], the results are compared to the LO band calculated by Platter et al. [||] 
and with our LO potential. The spread of the NLO values is not in conflict with EFT^ which 
predicts a 10% uncertainty at NLO. The upper (lower) boundary of the LO band was obtained by 
choosing different NN observables, a sA {B{d),a s ), to fit the LO LECs, and it was already pointed 
out in [Q| that this gives only a crude estimate of higher order effects. The difference between 
the 4 He binding energies obtained with the two potentials #2 and #6, which yield almost identical 
triton binding energies, serves as an estimate of the theoretical uncertainty. With this estimate and 
the points from the EFT^ potentials in Fig. [I] (left), we deduce a NLO correlation band of width of 
about 5 MeV centered around 28 MeV at the experimental B{t). The resultant prediction of 

B( 4 He) NL0 = (28 ± 2.5) MeV (4. 1) 

is consistent with the expected NLO uncertainty of about 10% and with experiment. The results of 
the AV18(+UIX) models lie within the proposed band as it is expected of all interaction models of 
at least NLO. 

Each of the nine squares in Fig. [I] corresponds to a set of A,'s with A^p = 0. Focusing on point 
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Figure 1: (left) Correlation between the triton and 4 He binding energies (Tjon line). The squares are results 
from the nine NLO EFT^ NN potentials. The solid (gray) black lines represent a smooth variation of A^p 
for a specific (LO) NLO potential. The upper (lower) dashed line was calculated at LO in ^ with LECs 
fitted to a Si t (B(d),a s ). The values for AV18(+UIX) are reported in JTofa, and those using a chiral expansion 
to NLO and N 2 LOin [§]. 

(right) The correlation between the triton binding energy and the real part of the spin singlet, S-wave scatter- 
ing length ao( 3 H e " n ) for elastic 3 He-n scattering. The dots represent the values of the NLO EFT^ potentials. 
The values for AV18(+UIX/IL2) [jTTJ] are results of RRGM calculations reported in [0. 



#1, the solid lines result from a smooth variation of A^p with the other LECs of the respective 
NLO and LO potentials held fixed. Both lines pass the datum well within the uncertainty bounds. 
The absolute value of A^p depends on the renormalization scheme chosen to obtain the two-body 
LECs. Our results also show at NLO that the two-nucleon force can be chosen such that A^p = 0. 
In general a non-zero three-nucleon counter term is mandatory for a correct renormalization of the 
triton. This term corrects an unphysical dependence of A > 2 low-energy observables on details of 
short-distance NN interactions, while low-energy NN predictions are unaffected. A three-nucleon 
interaction is therefore necessary to guarantee renormalization scheme independence of A > 2 low- 
energy observables. Once Aj^p is fitted to appropriate data, here B(t), all remaining dependence 
on the renormalization parameters, here A and input data, in B( 4 He) is a less than 10% higher 
order effect. This is pictured in Fig. [j] for two NLO and one LO potential, and we present it as 
evidence that no additional four-body parameter is needed at NLO, and hence the four-nucleon 
bound state is a universal consequence of the three-body system. Furthermore, the LO line (light 
gray), deviates more from the center of the correlation band, which is a confirmation of the order- 
by-order improvement of the accuracy of the EFT^ calculation. For a detailed analysis we refer to 
a forthcoming report [14]. 



4.2 Four nucleons: scattering 

In principle, all low-energy observables should be correlated with the triton binding energy 
via one parameter at NLO if the assertion that four-nucleon forces are at least N 2 LO is true. Also, 
the uncertainty in the prediction of other low-energy observables due to the truncation of the La- 
grangean should be comparable to the one identified in the previous subsection. We test both 
hypotheses by looking at the real part of the S-wave spin singlet scattering length ao( 3 He-n) for 
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elastic 3 He-neutron scattering. Its dependence on B(t) is shown for six NLO EFT^ potentials in- 
cluding the Coulomb interaction in the right panel of Fig. [I[ 

We find the naively expected decrease of Re{t?o( 3 He-n)} with increasing triton binding energy 
which maps out a band including the datum. Consistent with the analysis of the Tjon-line, the error 
estimate we deduce from the change in ao( 3 He-n) observed between the potentials #2 and #6 is 
about 10%. Again, the assumption of the center of the band at the experimental B(t) follows from 
the values in Fig. [I] (right). Thus, EFT^ reports 

Re{a ( 3 He-n)} NLO = (7.5±0.6) fm . (4.2) 

The NN model AVI 8 yields a value within the error band. Adding the UIX/IL2 three-body force 
(for details about those calculations see [ [l0| , [L^]) moves this point into the 10% NLO uncertainty 
range around the datum. This example of a low-energy scattering observable supports our previous 
conclusions: Every potential with the correct NN low-energy phase shifts and appropriately tuned 
three-body force, e.g., to give the correct triton binding energy, should not only predict the correct 
fi( 4 He) but also the experimental ao( 3 He-n) within a NLO error range of about 10%. No additional 
four-nucleon parameter is therefore necessary. 

4.3 Five and six nucleons 

The results presented in this subsection should give an outlook and serve as a demonstration 
of the feasibility of calculations in heavier systems with three-body- and Coulomb forces. Here, 
the 4 He-neutron ^ + )"P nase shifts were calculated as an intermediate step before looking at the 
6 He ground state. A bound state in this A = 5 channel would indicate a possible failure of EFT^ 
within its current ordering scheme. From the LO results shown in Fig. || (left) obtained with a 
cutoff A = 700 MeV and Ai^^nf fitted to a Sjt ,B(t), we conclude, preliminary, that no such state 
exists, and hence a modification of the power counting at LO is unnecessary. 
With this potential we calculated two observables in the 6 He system which we take as signatures of 
a halo structure. First, a relatively small charge radius (r 2 )*/ 2 compared to the matter radius (r 2 )l/ 2 , 
and second, a shallow bound state with respect to the 4 He-nn breakup threshold. The qualitative 
behavior of the two radii we calculated, 

(r 2 )i /2 = 3.261 fm , {r 2 )]' 2 = 5.678 fm , (4.3) 

hints at a halo structure. We attribute the quantitative discrepancy to data in part to the small model 
space used for the 4 He core in which it is only bound by fi( 4 He) 20 MeV. We expect a refinement 
of the variational basis to yield values consistent with EFT error margins. However, these margins, 
i. e. , the size of the expansion parameter of EFT^ in A = 6, have yet to be determined from the LO to 
NLO convergence and cannot be estimated a priori as the example of 4 He shows. In Fig. § (right) 
the dependence of (5( 4 He) - B( 6 He)) on B{t) is plotted for a varying A 3NF . For A = 700 MeV, the 
curve has a minimum close to the three-fermion system which corresponds to the physical triton. 
This means, that the halo structure is a peculiar property of nuclei, which is less pronounced in 
systems with the same two-body structure but different trimer properties. The halo character of the 
6-body state only follows for a certain range of the binding energy of the three-body bound state. 
It is not a universal consequence of the unnaturally large two-body scattering lengths alone. 
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Figure 2: (left) S-wave phase shift in the ^—channel for 4 He-n scattering. Compared are the results from 
a LO EFT^ calculation (red, solid line) to an R-matrix analysis (blue squares, [pj|). 

(right) Splitting between the 6 He and 4 He (lowest breakup threshold) binding energies. The thick (dashed) 
solid line was obtained with a LO EFTy potential with the two-nucleon force parameters held fixed while 
smoothly changing the three-body force. The set of RRGM parameters, the model space, was only optimized 
for the potential with A = 700 MeV. 



To confirm this finding, the aforementioned 4 He model space has to be enlarged to reach conver- 
gence in B( 4 He). A NLO calculation will yield an estimate of the expansion parameter of EFT^ in 
A = 6. At this stage of our analysis it is an open question whether or not 6 He is bound, a possible, 
relatively slow EFT convergence would also allow for an unbound nucleus. The value of the ex- 
pansion parameter will ultimately tell if the existence of a bound 6 He is a higher than NLO effect 
or not. Disregarding numerical issues related to the RRGM for the moment, the dependence for 
A = 800 MeV shows this possibility of an unbound system, i.e., with the change in the renormal- 
ization scheme, the shallow halo state (A = 700 MeV) disappears (A = 800 MeV). However, we 
attribute this observation primarily to the RRGM model space used, specifically the set of Gaus- 



sian parameters (see Eq. fl3.1|)) on the coordinates for the skin neutrons which were optimized for 
A = 700 MeV. From this analysis we deduce the following roadmap for ( 5 ' 6 ^He calculations to 
solidify our assertions. 

1. Check renormalization scheme independence by varying the cutoff (A G [300; 800] MeV) 
and by different input data (different NN P-waves, and A^p fitted to the charge radius of the 
triton instead of B(t) or the respective values of 4 He). 

2. For each potential thus fitted, the variational space has to be refined until full convergence is 
established. 



3. A NLO calculation with a clean, perturbative treatment of momentum dependent vertices 
has to be performed to demonstrate the convergence of the EFT expansion. 

4. Observables of subsystems have to be monitored carefully for consistency with the predicted 
theoretical uncertainty at the given order. For instance, the model spaces used for the 4 He 
cores do not meet this criterium and have to be refined as explained above. 
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5. Conclusions 

We presented an analysis of two universal correlations, the dependencies of the 4 He binding 
energy and a 3 He-neutron scattering length on the triton binding energy, with the nuclear Effective 
Field Theory "without pions". The calculations include a three-body force and the Coulomb 
interaction, and report the expected order-by-order improvement consistent with an expansion 
parameter p typ /A b « o. By that we demonstrated that a systematically improveable description of 
the 4 He-nucleus is possible with EFT^. In particular, we find no evidence that a four-body contact 
interaction is required to renormalize this system at next-to-leading order. The tetramer is a 
universal consequence of the trimer. 

We also presented preliminary results of an analysis in the five-nucleon scattering system, and of 
halo signatures of the 6 He system. The former reproduces the R-matrix data within leading order 
error margins while the latter yields the qualitative halo features but so far does not allow for final 
conclusions before certain numerical issues are carefully investigated. At this point of our 
analysis, the halo structure of 6 He is a universal consequence of two- and three-nucleon 
properties, while unnaturally large two-nucleon scattering lengths alone are not sufficient for a 
peculiar system like that. 

In a next step, we have to confirm the findings for the A > 4 systems at leading order by following 
the roadmap laid out in the previous section. Then, the theory will be used for predictions in 
nuclear systems with A < 8, including (un)stable Lithium and Beryllium states. We also aim 
towards an application to bosonic systems of atomic gases where recent measurements [|^] allow 
for an observation of the change of the renormalization scheme. As a prof of principle, a matching 



of EFT^ to a cluster effective field theory (e.g. [|15|]) would demonstrate this bottom-up approach 
of the construction of an EFT. 
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